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Abstract 

We consider the class ol d — A CFTs at finite temperature and chemi- 
cal potential that are holographically described within D — b Einstein- 
Maxwell theory with a Chern-Simons term. The high temperature phase, 
which is spatially homogeneous and isotropic, is dual to the AdS-Reissner- 
Nordstrom black brane solution. For sufficiently large Chern-Simons cou- 
pling, we construct new electrically charged AdS black hole solutions that 
are dual to the low temperature, spatially modulated phase. In this phase 
the current, associated with the abelian global symmetry, spontaneously 
acquires a helical order. The new black holes are stationary and also have 
Bianchi VIIq symmetry. 



1 Introduction 



Spatially modulated phases, in which the Euclidean spatial symmetry is sponta- 
neously broken down to some smaller subgroup, appear in condensed matter systems 
in a wide variety of settings, including spin density waves ^J, charge density waves [2] 
and FFLO states [3|[4], and are also anticipated in QCD at high baryonic density [s]. 
It is therefore of interest to investigate the properties of such phases for strongly 
coupled matter using the AdS/CFT correspondence. In this context it has been ar- 
gued that spatially modulated phases can arise for CFTs at finite temperature and 
charge density |6 11 and also in the presence of magnetic fields 12 ,13 . Other work, 
utilising the brane probe approximation, can be found in 14 17 . The simplicity of 
these constructions has led to the speculation that the typical ground states of holo- 
graphic matter at finite charge density and/or in a magnetic field could be spatially 
modulated (lOl. 

The first construction of black holes dual to spatially modulated phases was re- 
cently presented in the context of a D = 5 gravitational model with a gauge field and 



a charged two-form in 11 . Electrically charged AdSs black holes were constructed 
which are holographically dual to p-wave superconductors with a helical order. The 
black holes of |[llj are static and have a Bianchi VIIq symmetry, which is naturally 
associated with the helical order. It was also shown that at zero temperature the 
black holes become smooth domain wall solutions interpolating between AdS^ in the 
UV and a homogeneous but non isotropic ground state with a scaling symmetry in 



the IR, of a type similar to those found in 18 



Here we will consider another class of D = 5 models, first studied in this context 
in [g], namely Einstein-Maxwell theory with a Chern-Simons term. This class of 
models, parametrised by the strength of the Chern-Simons coupling, 7, can be used 
to study d = 4 CFTs with an abelian global symmetry, whose anomaly is fixed by 7. 
We will be interested in phases of the dual CFT at finite temperature and chemical 
potential with respect to the global symmetry which means that we need to construct 
electrically charged AdSs black holes. At high temperatures the CFTs are described 
by the standard AdS-Reissner-Nordstrom black brane solution. When 7 is greater 
than a specific critical value, 7c, the AdS-RN black brane has spatially modulated 
instabilities below a critical temperature, suggesting that the system moves into a 
spatially modulated phase in which the current acquires a helical order [6|. In the 
limit 7 — 00, where the back-reaction to gravity gets switched off, it was argued that 
the phase transition should be second order with mean field behaviour |7|. 
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The purpose of this paper is to construct the fully back reacted spatially modu- 



lated black holes for 7 > 7c. As in 11 , the key ingredient in the construction of the 
new black hole solutions is that they have a Bianchi VIIq symmetry associated with 
the helical order. With this symmetry, combined with time-translation invariance, 
we can use an ansatz involving several functions which just depend on a single radial 
coordinate. After substituting into the equations of motion we are led to a system 



of ODE's which we solve numerically. While the black holes in 11 were static here 
they are just stationary and this leads to the dual phase having spatially modulated 
momentum in addition to the spatially modulated pressure and shear that was seen 
in 11 . The new black hole solutions exist for temperatures lower than the critical 
temperature (which is not always the case |T^) and have less free energy than the 
AdS-RN black branch. The helical current phase is thus thermodynamically preferred 
and we show that it is always second order with mean field behaviour. We will see 
that the spatial modulation persists in the T — )■ limit and in this limit the entropy 
density approaches zero. 

Our black holes are dual to helical current phases which are reminiscent of the 
"chiral nematic (or "cholosteric" ) phase of liquid crystals (e.g. (20j). Recall that the 
order parameter for a nematic phase is a three-dimensional unit vector n, defined up 
to sign, called the "director" . In the chiral nematic phase there is a helical structure in 
which the director twists along an axis perpendicular to the direction of the director. 
For wave-number k the pitch of a helical phase is p = 2tt /k. While a general helical 
phase, including ours, is periodic with period p, for a chiral nematic it has period p/2 
because n = — n. One important issue is to calculate the temperature dependence of 
the pitch of the helical order (e.g. |2l]). In many materials the pitch increases with 
decreasing temperature but materials are known for which it decreases. We can obtain 
the precise temperature dependence of the pitch for our helical phases, finding that 
it monotonically increases, approaching a non-zero value at T = 0. Chiral nematics 
are well known to have interesting optical propertied such as selective reflection of 
circularly polarised light, and it will be interesting to explore analogues of them for 
our new black hole solutions using linear response theory. 



"'^This was recently explored using Lie algebra methods in 
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2 General setup 



We consider the D = 5 action given by 

S=jd'>x^ (i? + 12)-^F^,F^'^ -^JfaFAA, (2.1) 

where F = dA and 7 is a constant (the boundary terms will be given later). The 
corresponding equations of motion are given by 

If 1 _2 



d*F + ^FAF = 0. (2.2) 

These equations admit a unit radius AdS^ vacuum solution which is dual to a class 
of d = 4 CFTs with a global abelian symmetry. The metric, (7^1,, is dual to the d = 4 
energy momentum tensor, T™", and the gauge-field, A^, is dual to the d = 4 abelian 
current, J*". For example, when 7 = 2/\/3 ~ 1.1547 we obtain the bosonic content 
of D = 5 minimal gauged supergravity, and the class is known to include the most 



general class of = 1 SCFTs with type IIB or D = 11 supergravity duals 23 25 
We will be focussing on the range 7 > 7c, with 7^ ^ 1.1584. 

We will construct black hole solutions that are invariant under time translations 
and also have Bianchi VIIq symmetry. The Killing vectors associated with the latter 
are 8^2, d^^, which generate translations in the X2, X3 directions, respectively, and 
dxi — k{x2dx-^ — x^dx2)i where /c is a constant, which generates a helical motion con- 
sisting of a translation in the Xi direction combined with a simultaneous rotation in 
(x2,X3) plane. The corresponding invariant one- forms are given by 

uji = dxi, 

U2 = cos (kxi) dx2 — sin (kxi) dxs, 

UJ3 = sin (kxi) dx2 + cos (kxi) dx^, (2.3) 

which satisfy dui = 0, du2 = —kuji A and duj^ = kui A 002- The ansatz we shall 
consider is given by 

ds^ = -g f dt' + ^ + uj\ + r^e^" {uj2 + Qdtf + r^e'^'^ ujI , 

A = adt^huj2, (2.4) 

where /, (7, a, Q, a and h are functions of the radial coordinate r only. Note that 
when Q 7^ the spacetime is stationary but not static. The black hole event horizon. 
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located at r = where g{r^) = Q{r^) = a(r_|_) = 0, is, generically, the non-compact 
Lie group Bianchi VIIq. 

By substituting this ansatz into the equations of motion we find that / and g 
satisfy first order differential equations and that h, a, Q, a and b satisfy second order 
equations. Furthermore, these differential equations can be obtained from substitut- 



ing the ansatz directly into the action (2.1 ) and then varying the seven functions of r. 
The constant k is held fixed in these variations. As the expressions for the equations 
of motion are rather long, we just record the form of the action 

Q /"^s 2, J . ,f3f' 2h' r 
S = J d xr hf'i -g -g I — + + " 

_ [f'r^h + f [2rh + r^h') + / {r^h" + 2rh' + h)] 

, ,,2 2A;2sinh2(2a) e^^r''{Q'f Pr^e'^^g^ 
-2g[(x) p + + 

{p!f Qb'a' 1 fe-^"g 



+ 



+ 3 



2/2 p 2 V r2 




j d^xb{ba' -ab') . (2.5) 

It will be useful to observe that our ansatz, and hence the equations of motion, are 
left invariant under the following three scaling symmetries: 

r — )■ Ar, (t, X2, X3) — 7- A~"'^(t, a;2, xs), g ^ X^g, a — )■ Aa, b ^ Xb; 
xi — )■ A^^xi, h -> A/i, k ^ Xk; 

t^Xt, /^A-V, a^A-^a, Q^X-^Q; (2.6) 

where A is a constant. 

The equations of motion admit the electrically charged AdS-Reissner-Nordstrom 
black brane solution. It has h = r, f = l,a = Q = b = 0, and hence, 

ds^ = -gdt^ H h [dxj + dxj + dxf) , A = adt, (2.7) 



with 



The AdS-RN black brane is static and has Euclidean, IS0{3), symmetry. It has 
temperature T = (6r^ — fi^)/67Tr^ and describes the high temperature, spatially 
homogeneous and isotropic phase of the dual CFTs when held at finite chemical 
potential with respect to the global abelian symmetry. 
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2.1 Asymptotic AdS^ and near- horizon expansions 

We will be interested in new black hole solutions that asymptotically approach AdS^ 
in the UV and are dual to d = 4 phases where the breaking of the Euclidean symmetry 
to a helical order is spontaneously generated. By analysing the equations of motion 
we can construct the following asymptotic expansion as r — cxd: 

^4 ' 



/ = /o (] 



h = r (1 + % + 



^4 



<3 = /o(^ + -- 

a = fo (/^ + ^ + 



^ + --- (2-9) 



At a convenient juncture we will use the scaling symmetries (2.6) to set /o = yU 



The UV data is then specified by seven parameters M,Ch,Ca,CQ,q,Cb and k. Note 



that we have fixed the asymptotic fall-off of h in (2.9) so we can no longer use (2.6) 
to scale k. The holographic interpretation of these parameters will be discussed later. 

At the black hole horizon, located at r = r_|_, the functions have the analytic 
expansion 

9 = 9+ (r - r+) + . . . 

/ = /+ + ••• 
h = h+ + ... 
a = a+ + . . . 

Q = Q+ir-r+) + ... 
a = (r — r+) + . . . 

b = b+ + ... (2.10) 

Regularity of the metric at the black hole horizon can easily be seen by using the 
in-going Eddington-Finkelstein coordinates v,r where v ^ t + {g^f^)^^ ln{r — r+). 
The full expansion is fixed in terms of the seven constants /+, «+, h^, Q^, a_|_, 
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and r_|_. In particular, the coefficient g+ is fixed by these constants: 

+ 4-^ r+ (2.11) 



After fixing the scahng symmetries (2.6) we have seven UV parameters and seven 



IR parameters. We have two first order differential equations and five which are 
second order, so a solution is fixed by twelve parameters. Thus, generically, we 
expect a two parameter family of black hole solutions, which we will label by k and 
temperature T. 

2.2 Thermodynamics 

To analyse the thermodynamics of the black hole solutions we will need to calculate 
the on-shell Euclidean action. Additional details are presented in appendix |A} We 
analytically continue by setting t = —it and in order to get a real metric and vector 
field we should also write (5+ = cl+ = id^. Near r = r+ the Euclidean solution 
then takes the approximate form 



dsl ^ flir - r^) (dr + ^2) + 

\ 9+J+ J 



+ hi ul + r2 e2"+ [u^f + r2 e-2-+ ujI , 
AK,a+{r-r+)dT + buj2. (2.12) 

Regularity of the solution at r = r+ is then easily seen by making the coordinate 
change p = 2g^^^'^{r — r+)^/^ and making r periodic with period At = 47r/((yf+/+), 
corresponding to temperature T = (/qAt)^^. We can also read off the area of the 
event horizon and, since we are working in units with IGttG = 1, we deduce that the 
entropy density is given by 

s = 4:7ir^h+. (2.13) 
We will consider the total Euclidean action, Ixot, defined as 

iTot = I + hndy , (2.14) 



where / = —iS and hndy is the Euclidean boundary action of 26 , including counter- 
terms, which is given explicitly in appendix |Aj We next define the potential W, and 
a corresponding density w, for the grand canonical ensemble via W = T[lTot]os = 
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wvol^, where [Irotlos is the on-shell Euchdean action and V0I3 = J dx^dx'^dx^. A 
calculation reveals that w can be expressed in two equivalent ways 



w 



-M = 3M + 8ch + 2fxq - Ts 



(2.15) 



with the equality of the two expressions giving a Smarr type formula. 

A variation of the bulk action / gives equations of motion and boundary terms. 
Thus an on-shell variation only gets contributions from the boundary. We hold k 
fixed in these variations, for reasons we discuss in the next subsection, and for the 
Euclidean black hole we then only get contributions at r — )■ 00. Combining this with 
an on-shell variation of the boundary action Ihndy and using the asymptotic expansion 



(2.9) we deduce that w = w{T,fi) and 

6w = —s6T + 2q6fi . 



(2.16) 



To illuminate the holographic meaning of the constants appearing in the UV 



expansion (2.9), we now compute the expectation value of boundary stress-energy 



tensor and the current. The relevant terms for the stress tensor are given by [27] 



lim r' 



[-2Kmn + 2{K - 3)(^oo)mn + 



(2.17) 



Using (2.9) we obtain 



{Tu)= fl (3M + 8C,) 
(Ttej) = 4/oCq cos(A;a;i) 
(Tteg) = -4:foCQ sin(A;xi) 
(T,,,,) = {M + 8ch) 
(T,,,,) = (M + 8c„ cos(2A;xi)) 
{Txsxa) = (M - 8ca COS {2kxi)) 
{Txixa) = -8ca sin(2/i;xi) . 



(2.18) 



One can check that this is traceless {goo)"^"" {Tmn) = 0. Setting /q = 1 we see that the 
energy density of our solutions, e, is given by 



e = 3M + 8ch 



(2.19) 



Furthermore, we also deduce from (2.15) and (2.16) that the first law can also be 
written in the form: 



6e = T5s - 2fi5q . 



(2.20) 
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Returning to (2.18) we see that Cq, specifies spatially modulated pressures and shear in 
the {x2, X3) plane, with the length scale of the modulation fixed by the wave-number 
k. The pressure in the Xi direction is given by M + Sch- If we define p to be the 
average of the three pressures we have p = M + 8/3c/j, and the Smarr formulae in 



(2.15) can be written e + 2nq = Ts + p — 8/3ch- These features were also seen for the 



black holes found in 11 . A new feature of the black holes we are considering here 
is that there is spatially modulated momentum in the {x2,X3) plane specified by cq 
and k. 

We next calculate the expectation value of the current. The relevant terms are 

(Jm) = lim r'[Frm + ...]. (2.21) 



given by 26 



where the ellipsis refer to terms that will not be relevant here. Using (2.9) we obtain 



(Jt) = -2/og 
(^xi) = 

(Jxi) = -2cfeCos(/cxi) 

{J^,) = 2ctsm{kxi). (2.22) 

From the temporal component we see that the constant q fixes the charge density. 
From the spatial components we see that fixes the strength of the spontaneously 
generated spatially modulated helical current. It is clearly circularly polarised. 

2.3 Variations of k 



In the variations to get (2.16), or equivalently (2.20), we held the wave-number k 
fixed. One can consider arbitrary non-normalisable and normalisable deformations of 
the fields and then expand them in a complete basis of functions. Here we are viewing 
k as labelling one of the modes and hence should not be varied to obtain the equations 
of motioi]|^ Furthermore, it should not be varied to obtain an on-shell variation of 
the action in order to obtain the thermodynamics for a particular solution labelled by 
k. As we discuss further in the next section, we will obtain a two parameter family of 



^An analogous procedure was employed in [ll] and also, essentially in a field theory context, 
in I?]. To clarify this point, in appendix |A] we consider a more general set-up in which we also allow 
a more general deformation parameter at infinity. Specifically, we consider 6 — /i;, + ^ + . . ., with 
fib and Cb being the non-normalisable and normalisable deformations of the magnetic part of the 
gauge field for a particular mode k. 
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black hole solutions to the equations of motion that depend on k and T (see figure |2|. 
At fixed temperature, this should be viewed as a moduli space of solutions, labelled 
by k and we should choose the solution labelled by a specific value of k that has the 
smallest free energy w. This leads to a one-parameter family of thermodynamically 
preferred solutions labelled by T, given by the red line in figure |2] In fact this red 
line is specified by the condition that the action Ixot is stationary with respect to a 
free variation of k: 

dr{k(^ h + — J-3^(^"-"^)^=°- ^'•23) 



In appendix |A] we will discuss how this arises from contributions to varying the action 
at xi = ±00. 

Another perspective is to consider the xi direction to be periodic with xi = Xi + L. 
In this case we should only consider discrete wave- numbers k = n/{2TTL), for arbitrary 
integer n, and there is no issue of varying k to obtain the equations of motion. In this 
case there will be a discrete set of solutions on figure |2] and, at a fixed temperature, 
one should just choose the one with smallest free energy, as usual. One finds that 
the system will, in general, jump discontinuously from one branch to another giving 
a series of first order phase transitions. In the limit that L — )■ 00 we will recover the 
continuum picture that we have discussed above. 

It is worth noting that the same kinds of issues also arise for homogeneous and 
isotropic phases. For example, recall the basic s-wave holographic superconducting 



black holes 28 , 29 . In this setting the AdS-RN black brane, which describes the 
high temperature phase, becomes unstable to the formation of charged scalar hair. 
Although only scalar modes with k = have been discussed in the literature, a 
linearised analysis for k will produce a curve analogous to that in figure [T] but 
it will now be symmetric about k = 0. Hence, there should also be a two-parameter 
family of superconducting black hole solutions labelled by T and k. In this case, 
however, the thermodynamically preferred curve of solutions will just be the solution 
with k = 0. 



3 Helical black holes 

Based on the analysis of linearised perturbations about the AdS-RN black brane 
solution carried out in [6] we expect to be able to construct spatially modulated black 
hole solutions provided that the Chern-Simons coupling 7 is larger than 7^ ~ 1.1584. 



9 



Figure 1: The curve denotes the critical temperature at which the AdS-RN black 
brane becomes unstable and also where the new branches of helical black holes, given 
in figure [2| appear. The plot is for 7 = 1.7 and ^ = 1. 



We will now set 7 = 1.7, and hence 7/7C ~ 1.47, but we have checked that several 
other values lead to qualitatively similar results. For this value the linearised analysis 
of [g], which we summarise in appendix [B| leads to the curve presented in figure [l] 
which denotes, for a given value of k, the temperature at which the AdS-RN black 
brane becomes unstable. Hence for k in the range 0.47 ^ k < 3.05 we expect to be 
able to find the new black hole solutions. 

The new helical black hole solutions are obtained by solving the equations of 



motion numerically for the ansatz (2.4) with boundary conditions at the asymptotic 



AdS^ boundary given in (2.9), and at the black hole horizon given in (2.10). We 



use the scaling symmetries (2.6) to set fo = fx = 1. As mentioned earlier a simple 
parameter count indicates that we expect, generically, a two-parameter family of 
solutions which we take to be labelled by temperature T and wave-number k. In 
practise we fix a specific value of k and then construct a one parameter family of 
solutions labelled by the temperature T. We considered twenty different values of k, 
in the range 0.6 < A; < 1.8 (focussing on the peak of the curve in figure [T]), and we 
have displayed our results in figures [2] - 14} 

Figure [2] shows the two-parameter family of solutions and their free energy w. We 
first note that the boundary of the surface projected onto the {k, T) plane reproduces 
the curve of critical temperatures as a function of k where the AdS-RN black brane 
becomes unstable given in figure [l] We next note that for any fixed temperature the 
helical black holes have less free energy than the AdS-RN black hole for any value 
of k. Thus, from figure |2] we deduce that there is a second order phase transition at 
T = Tc ^ 0.0627 at A; = fcc ~ 1-32 with the system moving from a homogeneous and 
isotropic phase to a spatially modulated helical phase. As the temperature is lowered 
we need to find the value of k for which the black hole has the lowest free energy. This 
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Figure 2: The two-parameter family of helical black holes, labelled by temperature 
T and wave- number and their free-energy w. The red line denotes the thermody- 
namically preferred locus, which minimises w over the moduli space of solutions at 
fixed T labelled by k. The plot is for 7 = 1.7 and n = 1. 




Figure 3: Plot of Ch versus k for the one parameter family of helical black hole solu- 
tions given in figure [2] for the representative temperature T ^ 0.0535. In particular 
Ch = on the thermodynamically preferred red line of solutions given in [2] The plot 
is for 7 = 1.7 and fi = 1. 



leads to the one-parameter family of thermodynamically preferred black hole solutions 
which are marked with a red line in figure 1. We have also checked numerically that 



this red line coincides with imposing (2.23). Interestingly, our numerical analysis 
indicates that all of the black hole solutions on the red line have c/j = (a similar 
phenomenon was also seen in [ll]). In figure [3] we have plotted the behaviour of Ch 
versus k for a representative temperature of T 0.535 and we see that it vanishes 
along the red curve as well as on on the boundary curve of figure |2} It would be 
interesting to understand the underlying reason for this behaviour. 

In figure |4] we show the behaviour of various other physical quantities as a function 
of temperature for the thermodynamically preferred branch, given by the red line in 
figure [2| and marked with red lines in figure |4| as well as for the AdS-RN black hole 
solution, marked with blue lines in figure |4} The first two panels show the free energy 
w for both solutions and the wave-number k for the red line. The pitch, p of the 
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Figure 4: The red lines plot various physical quantities against temperature T for 
the thermodynamically preferred helical black hole solutions on the red line in figure 
|2} The blue line refers to the AdS-RN black hole solution, w is the free energy and 
k is the wave-number of the helical order, cq and Cq fix the spatially modulated 
momentum and stress/strain in the {x2iX^) plane, respectively, q and Cb determine 
the size of the charge and the spatially modulated current, respectively, and s is the 
entropy density. The plots are for 7 = 1.7 and fi = 1. 



helical order is given by p = 271 /k and hence we see that the pitch monotonically 
increases as the temperature is decreased. The next two panels show the behaviour 



of Cq and Ca which, we recall from (2.18), determine the strength of the spatially 
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modulated momentum and stress/strain in the {x2,X3) plane, respectively. The next 
two panels show the behaviour of the charge density q and also the behaviour of 



Cfc, which we recall from (2.22), determines the strength of the spatially modulated 
helical current. The final panel shows the behaviour of the entropy density, s. By 
analysing the behaviour close to the critical temperature, T = T^. ^ 0.0627, we find 
the following mean-field behaviour 

CQ^ 2953T,^ , c„^107T,^(^l-| 

Cfe^SaOT^ A: ^T, (^11.5 + 9.58^^ . (3.1) 

3.1 Low temperature behaviour 

It is clear from figures |2] and |4] that the helical order persists as T — )■ 0, and in 
particular k approaches a non-zero value in this limit. Furthermore, the entropy 
density approaches zero. However, we have not yet been able to pin down the precise 
behaviour of the solution in this limit, and we leave this interesting issue for further 
instigation. However, we record a couple of conclusions based on our numerical 



results. Amongst the coefficients in the near horizon expansion (2.10), we find that 
/i+ — > oo and a+ — t- — oo, with the entropy density, s = inr'^h^, going to zero and 
r+/i+e~"+ approaching a constant value. We also find that /+, a+ vanish while 
go to constant values. Starting from the phase transition the value of F^^^F'^'^ at the 
horizon is — and this monotonically increases and approaches 24 as T — )■ 0. Thus 
associated with the vanishing entropy density, we have an expulsion of electric charge 
somewhat reminiscent of [30]. We also evaluated the horizon value of the Ricci scalar 
and the square of the Ricci tensor and it appears that i? — t- — 18 and R^^R^'^ 108. 



It is tantalising that that these are the same values for the Schrodinger solution of 31 . 



4 Final Comments 

We have constructed, numerically, a new class of electrically charged AdSs black 
holes that are dual to = 4 CFTs at finite charge density acquiring a helical current 
order via a second order phase transition. We have extracted a number of physical 
properties of the helical phase including the temperature dependence of the wave- 
number k that fixes the pitch, p = 27: /k, of the helix. We have shown that the 
pitch monotonically increases as the temperature is lowered but approaches a finite 
value as T — > 0. Furthermore, our numerical results indicate that the entropy density 
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goes to zero in this limit. It will be very interesting to further analyse the precise 
behaviour of our solutions as T — )■ in order to better understand the emergent 
spatially modulated ground state at T = 0. 

Another direction is to calculate various transport coefficients by calculating var- 
ious two point functions. There will be a number of different channels to analyse 
and we expect a rich structure. It will also be interesting to explore the related 
hydrodynamics of the black holes. 

Our numerical results imply that the thermodynamically preferred helical black 
holes (the red line in figure [2]) have the property that the expansion coefficient Ch, 



appearing in (2.9) and entering the definition of the energy density and the pressure 



in the direction of the axis of the helix (see (2.18)), is exactly zero. The same 
phenomenon was also seen for the helical superconducting black holes of |11!| and it 
is desirable to have a better understanding of this property. 

The construction of the D = 5 black hole solutions describing spatially modulated 



phases here and in 11 has been facilitated by the fact that they are static and 
have a Bianchi VIIq symmetry. This leads to the construction of a co-homogeneity 
one ansatz for the D = 5 fields and hence solving ordinary differential equations. 
Moving to D > 5 (obviously of less interest in making connections with condensed 
matter systems), many generalisations are possible while staying within the realm of 
solving ODEs. However, constructing D = 4 black holes that are dual to spatially 
modulated phases (as in the models [8|[T2]) will necessarily involve solving partial 
differential equations. Wile this is technically more challenging we expect to see new 
phenomenon. 
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A More on thermodynamics 



Here we will expand upon the discussion of the thermodynamics that we summarised 
in section 2.2[ We found it illuminating to consider this issue in a slightly more 
general setting in which we allow for a non-normalisable fall-off in the magnetic part 
of the gauge field at infinity: 

6 = /ife + ^ + .... (A.l) 

We emphasise that our new helical black hole solutions all have /ib = 0. The asymp- 
totic expansion at infinity now reads 

2 (. _M k^^ij + 16cg + _ 2A;VbCblnr fcVg Inr^ 

f=f(^. + ^ _ fcVfcln^ , -49fcVg + 48fcVfcCfc-432c,^ 
J J^y+ ^4 i6r4 ^ 1728r6 

(-/tVfc + 18A;V6C6)lnr _ A;Vb(lnr)2 

/ Ch k'^nllnr 35A;Vb ~ 96fc^/ifeC6 144c^ {k'^fi^ — 3k'^fibCb)lT^r 
V ^ ^ ^ 16r4 + 1728^^6 + ^6 

fcV^(lnr)^ 
48^6 

_ c« A; Vft In r 576cak^ - 59k^fil + 96k'^fibCb - UAcl 
" ~ ^ 16r4 1728r6 

(-7A:Vb + Inr _ A;Vfe(lnr)2 

+ 144r6 48r6 " ' 

n - f -I- ^Vfeg ~ ^'^^bq + Sk'^CQ k'^fibq Inr 

, g --fk^^l + 8'ykfibCb 'jk^fillnr 
^-fo{f^ + -2+ 8;:^ + --- 



Cb k'^fiblnr -3k'^Hb + 8k'^Cb + l6'yknbq /cVfelnr , _^ ^, 



The expansion at the black hole horizon is presented in (2.10). 
We will consider the total Euclidean action, Ixot, defined as 

ItoI = I + hndy , (A- 3) 

where / = —iS and the (standard) Euclidean boundary action, hndy, is given by an 
integral on the boundary r — )■ oo [26] : 

hndy = j dr(fx |^_2i^ + 6 - ^ In r F^^F"^" + . . . (A.4) 
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Here K = g'^'^Vmf^n is the trace of the extrinsic curvature of the boundary, where 
is an outward pointing normal vector, and Qoo is the determinant of the induced 
metric. The Inr term is required to remove the divergence associated with the trace 
anomaly T^™ = —^FmnF'^^ and the ellipsis refers to a Ricci scalar term which will 
not be relevant for the ansatz and boundary conditions that we are considering. For 
our ansatz we have 



' hndy 



volaAr lim r hfg 

r— >oo 



1/2 



Inr 

2 



a' - Qb'^^ 



P 



9{b' 



(A.5) 



where V0I3 = J dx^dx^dx^ . We next point out two equivalent ways to write the bulk 
part of our Euclidean action on-shell: 



I OS = volsAr 



dr 



r+ 

00 



2rghf + -^pQQ' + ^he-'Vgbb' + —r'h {a' - Qb') bQ + -kjab^ 



vol::iAr 



dr 



r+ 



r^hfg' + 2r^hgf - -r^e^'^QQ' - K^ha {a' - Qb') - h^ab^ 
J J 3 

(A.6) 

Notice that the first expression only receives contributions from the boundary at 
r — )■ cxD since g{r^) = Q{r^) = a(r+) = 0, while the second expression also receives 



contributions from r = r_|_. Using the expansions at the AdS boundary (A. 2) and at 

the black hole horizon (2.10), and combining with (A.5) we obtain the following two 

equivalent expressions for the total on-shell action: 

1 11 
[Ixoths = vol3-[-M - fibCb- —d^^ + -fJ'fJ'lk-f] 
1 12 o 



Volq 



'T 



1 



3M + d>Ch + 2/ig -Ts- -^Ik"^ - -\i\iik-i 



(A.7) 



A variation of the bulk action / gives equations of motion and boundary terms. 
Thus an on-shell variation only gets contributions from the boundary. We hold k 
fixed in these variations and then we only get contributions at r — )■ 00. Combining 
this with an on-shell variation of the boundary action Ihndy and using the asymptotic 
expansion (A. 2) we eventually obtain 

{^Ito^os =vol3Ar ( 8c;, + 3M + 2/ig 



1 .2,2 1 



[i^k - -///ift/cT 5/0 + /o 2g - -/ifc/cT b[i 



1 



(A.8) 
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In this variation we are holding Ar fixed and hence AtSJq = —T~'^ST. We next define 
the potential W, and a corresponding density w, for the grand canonical ensemble 
via W = T[lTot]os = WV0I3. We deduce that w = w(T,fi,fii,) with 



Sw 



-s6T+ [ 2q - -nlk-f ) 6fi 



1 ,2 1 , 

2c6 + -yUfefc - -fifXbk'j 



(A.9) 



We now compute the expectation value of the stress tensor and the current. For 



the former we have 26 



(T^n) = lim 



~l" 2(7^^ 3)((7oo)mn ~l~ ( F-rn' Fnp ^i,9oo)'mnFpqF^'^ J In T -|- . . 

(A.IO) 



Using our expansion at the AdS boundary (A. 2) we obtain 

{Tu)= fl (^3M + 8ch-^filk' 

(Ttxi) = 4/oCq cos(A;xi) 
(^^3) = -4/oCq sin(fca;i) 

(T,,,,) = M+ fsc, + ^ filk^] cos i2kxi) 

1 



(T,3,3) = M - (8ca + - fitk'j cos (2/cxi) 
(^^2x3) = - ( + ^fJ'lkA sin(2fca;i) . 



(A.ll) 



Observe that {Tm^) = —filk /Q = —j^FmnF'^"', as expected (here we are raising 



indices with and the factor appears because it also appears in (A.IO)) . Setting 
/o = 1 we see that the energy density is given by 

1 



2;.2 



e = 3M + 8ch- - /ifefc 



(A. 12) 



Observe that the equality of the two expression in (|A.7|) imply the Smarr type formula 

(A. 13) 



If we define the average pressure p = {{Tr^-^x^) + {Txixi) + (^xixi))/3, this can also be 
written in the form 



e + p = sT - 2fiq + fc^ - fibCh - ^f4k'^ + ll^fAf^l 



(A. 14) 
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We next calculate the expectation value of the current. The relevant terms are given 
by [26] 



lim 

r—>OD 



r^F '^^'^A F + V F" lnr + 

' rm Q I "^m ^n-' pq ~ ^ n-'- m ' ~ ■ ■ ■ 



(A. 15) 



where V is the Levi-Civita covariant derivative with respect to the boundary metric 



goo- Using the expansion (A. 2) we obtain 

1 



(Jt) = - 
(J.,) = 



{Jxi) = - ( 2c6 + ^/ifc/i;^ - ^/i/i6^7 ) cos(/cxi) 



{Jxs) 



'^Cb + ^fJ^bk^ - ^fJ'fJ'bkj ] sin(A;xi) . 



(A.16) 



In terms of the current, the Smarr formula (A. 14) can be written in the form 



(A.17) 



e + p = sT + fi{Jt) + |c/i + {{Jxi) cos /cxi - (J^^g) sin/cxi) + ^^ilk"^ 

As we discussed in section 2.3 we view k as labelling a particular mode and hence 
should not be varied. To amplify this point, it is useful to refer back to our ansatz 



(2.4) and define the xi dependent variation [5y4(fc)]m via 

[5A{k)\ = 5^^ 

mk)u = 

[6A{k)]^^ = 6fib cos(/i;xi) 
[6A{k)],^^ = -6fib sin(fca;i) . 



We then find that the first law (A. 9) can be written in the form 



SW^Jdx.dx.dx.{-sST^iJ-i-kmSAik)^^), 



(A. 18) 



(A.19) 



where (J'"(— A;)) is given in (A.16) and we note the integrand is actually independent 



of Xi. In particular, we interpret 6 fib as parametrising a specific mode, labelled by k, 
of a non-normalisable deformation of the gauge-field, 6A{k). 



A.l Another perspective on (2.23) 



Let us consider a variation of the gauge field part of the bulk action (2.1 ). This gives 
the boundary term 



6S, 



gauge 



dM 



*F + -AAF 
3 



A6A 



(A.20) 
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Here we would like to focus on variations of wave-number k which give rise to a 
contribution at the boundary of the non-compact Xi direction. Given our ansatz 



(2.4) we have 



6A 



-bxi 003 6k 



(A.21) 



We take the boundary to be at xi = ±L/2 and then take L — )■ 00. A calculation 
shows that the relevant part of the integrand is 



*f + -aaf 

3 



A5A 



-dt A dxo A dx'i A dr 



h 



— '^b{ba' — ab')j Xi5k + . . . 

(A.22) 



After evaluating this at Xi = ±L/2 and then dividing by L in order to find the 



variation of the density, we are led to the third and fourth terms of (2.23) (taking 
into account a minus sign since here we are looking at the Minkowski signature space- 
time). 

If we now vary the Einstein-Hilbert term we get the boundary term 



6S 



EH 



JdM 



(A.23) 



where 7 is the induced metric on the boundary dM and for the boundary compo- 
nents defined by xi = ±L/2 we have that the unit normal vector is n = dx^- 



Substituting a variation of the metric obtained by varying k in our ansatz (2.4), we 
obtain contributions at xi = ±00 which lead to the first and second terms in our 



formula (2.23), again up to a minus sign. 



This calculation show that variations of k in our ansatz (2.4) are associated with 
boundary contributions to the variation of the action at Xi = ±00. Since we do not 
want to modify boundary conditions at Xi = ±00, k is a parameter to be held fixed 
in obtaining the relevant equations of motion. 



B Linearised Analysis 

We summarise the analysis of linearised perturbations about the AdS-RN black brane 
solution considered in [6] in the language of this paper. Specifically, we consider the 
perturbation 

Q^eQ, b^eb (B.l) 
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with h 



a 



in (2.4) around the AdS-RN black brane solution (2.7), for small e. 



At first order in e we obtain two coupled ODEs, linear in Q and b, given by 

r-^ [r^Q')' - k^r~^g-^ Q + r'^ a'b' = , 
r~^9~^ {rgb')' - k\-^g-^ b + r^g-^a' Q' + kr-^g-^-fa' 6 = 0, (B.2) 

with g and a as given in (2.8). To make contact with equation (4.17) of f6] one should 



make the following identifications: u = r^/r, i}j{u) = —\/?)r^Q'{r)/rj^^ b{r) = (p{u), 
q = ii/{r+y/3), kthere = khere/r+, a = 'J / A, q = /i/(r+A/3). 

In order to find a normalisable static linearised mode of interest we should impose 
the following boundary conditions. At the horizon, r — )■ r_|_, we demand that 



Q = Q(+) (r - r+) + . . . 
b = b+ + ... . 



(B.3) 



We are only interested in deformations of the CFT given by the temperature T and 
the chemical potential fi. Hence as r — t- oo we demand that 



b = — + 



[B.4) 



A solution to B.2 is specified by four integration constants and hence, for a given T, fi, 



we expect a unique solution (if any). By numerically solving (B.2) we find solutions 



that are summarised in figure [T| for the special value 7 = 1.7 (and also /i = 1), which 
agrees well with figure 2 of (6). 
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